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Abstract 

The two-dimensional superconductor-insulator transition in disordered ultrathin amorphous bismuth films has been tuned 
both by electrostatic electron doping using the electric field effect and by the application of parallel magnetic fields. Electrostatic 
doping was carried out in both zero and nonzero magnetic fields, and magnetic tuning was conducted at multiple strengths of 
electrostatically induced superconductivity. The various transitions were analyzed using finite size scaling to determine their 
critical exponent products. For the electrostatically tuned transition the exponent product vz — 0.7 ±0.1, using data from 
intermediate temperatures down to 60 mK. Here v is the correlation length exponent and z is the dynamical critical exponent. In 
the case of electrostatically tuned transitions in field, and the field-tuned transtions at various values of electrostatically induced 
superconductivity, from intermediate temperatures down to about 100 to 150 mK scaling was successful with vz — 0.65 ±0.1. 
The parallel critical magnetic field, B c , increased with electron transfer as (An c — An) 0-33 , and the critical resistance decreased 
linearly with An. However at lower temperatures, in the insulating regime, the resistance became larger than expected from 
extrapolation of its temperature dependence at higher temperatures, and scaling failed. These observations imply that although 
the electrostatic- and parallel magnetic field- tuned superconductor-insulator transitions would appear to belong to the same 
universality class and to be delineated by a robust phase boundary that can be crossed either by tuning An or B, in the case 
of the field-tuned transition at the lowest temperatures, some different type of physical behavior turns on in the insulating 
regime. 



I. INTRODUCTION 

Continuous quantum phase transitions are transitions 
at absolute zero in which the ground state of a system is 
changed by varying a parameter of the Hamiltoniani^ 
The transitions between superconducting and insulat- 
ing behavior in two-dimensional superconductors tuned 
by magnetic field or disorder (thickness) are believed 
to be such transitions. Early experiments and theo- 
ries seemed to support a picture of only two ground 
states. Historically the first theoretical approach to 
the superconductor-insulator (SI) transition was based 
on Cooper pairing being suppressed by the enhancement 
of the Coulomb repulsion between electrons with increas- 
ing disorder ^£ In effect the order parameter would be 
suppressed to zero in the insulating regime. The absence 
of a gap in the density of states in tunneling studies of the 
insulating regime has been interpreted as evidence for a 
zero superconducting order parameter amplitude in the 
insulating regime. 6 Another approach to the SI transi- 
tion was based on the transition being governed by phase 
fluctuations In this "dirty Boson model," the insu- 
lator is a vortex condensate with localized Cooper pairs, 
in contrast with the superconductor, which is a Cooper 
pair condensate with localized (pinned) vortices. Elabo- 
rations on this model have included effects of fermionic 
degrees of freedom. 10 If phase fluctuations are pair break- 
ing and destroy the superconducting energy gap, then 
this picture would also be consistent with the tunneling 
studies. 

The three conventional approaches to the tuning of SI 
transitions all involve uncontrolled aspects of morphol- 
ogy. In the first method, a relatively thick (^ 100A ) 
superconducting film is produced and magnetic fields are 
applied to quench the superconductivity. These films are 



two dimensional in the sense that the coherence length 
and penetration depth are much greater than the films' 
thicknesses. In this approach, the strength of the vor- 
tex pinning, which depends on the nature of the disor- 
der, may determine the outcome of the measurements. 
In some measurements on films with relatively low sheet 
resistances^, the transition can also be described using 
quantum corrections to conductivity^ and may not be 
quantum critical at all. 13 In a second approach, the thick- 
ness of a film is increased in small increments, tuning 
from insulator to superconductor However, films of 
different thicknesses may have different morphologies. In 
a third method, which has been used in some studies 
of In2 03 films, an insulating film is thermally annealed 
to produce superconductivity^*^ However, thermal an- 
nealing may involve the alteration of morphology and 
chemical composition. 

Successful finite-size scaling analyses with film thick- 
ness or magnetic field as tuning parameters have resulted 
in critical exponent products, vz, in the range of 1.2 to 
1.4^A*A2*A2i2fii2i which have been suggested to result from 
the transition being dominated by percolative effect s^2iSi 
as this number is close to the exponent in 2D percolation. 
One investigation of a perpendicular field tuned transi- 
tion has yielded 0.7 as the product. 19 

The simple two-ground state picture has been chal- 
lenged in recent work, which appears to indicate that 
there is an extended intermediate metallic regime over 
an extended range of the tuning parameter. The physi- 
cal evidence for such a regime is that resistances become 
independent of temperature at the lowest temperatures 
both in perpendicular field tuned transitions 24 and in 
thickness-tuned transitions. 25 In the case of the magnetic 
field experiments, the metallic regime abruptly disap- 
pears as the magnetic field is reduced, leading to true 
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superconductivity. This is further evidence of possible 
sample inhomogeneity. There are several theories that 
describe an intermediate metallic regime2£i 2 i2& in homo- 
geneous samples. However it is not established that this 
regime is indeed intrinsic, and is not a consequence of ei- 
ther sample inhomogeneity, or the failure to cool a film. 
In overtly granular systems there is a clear intermediate 
metallic regime that is found at temperatures in excess 
of 1 K that is probably intrinsic 

Another phenomenon reported by a number of groups 
is the appearance of a large peak in the resistance at 
fields in excess of the critical field for the SI tran- 
sition for superconducting amorphous 1^03 and TiN 
films j 13 i 30 i 31 i 32 i 33 i 34 i 35 i 36 or intrinsically insulating amor- 
phous Be films. 37 With the exception of the work of Gant- 
makher et a/., 38 most of these studies have been con- 
ducted in perpendicular magnetic fields No quantitative 
theoretical explanation has yet been put forward to ex- 
plain these very large values of resistance found on the in- 
sulating side of the SI transition. It is not clear that this 
phenomenon is an intrinsic property of a homogeneous 
material or results from some mesoscale inhomogeneity. 

Since there are important concerns about morphology 
and disorder in the above-mentioned studies, the aim of 
the present work was to attempt to clarify these issues 
by studying SI transitions in which the level of physi- 
cal disorder was fixed, and in which the outcome was 
not dependent upon the degree of vortex pinning. This 
is possible by inducing superconductivity in an insulator 
by electrostatic doping using the electric field effect 39 and 
then applying a parallel magnetic field to drive the film 
back into the insulating state. The same level of chemical 
and physical disorder may be shared by both the intrin- 
sic insulating state, and the insulating state in which the 
electrostatically induced superconductivity is quenched 
by magnetic field. We will present arguments to the effect 
that electrostatic doping does not alter physical or chem- 
ical disorder, but changes the coupling constant that de- 
termines superconductivity. Parallel magnetic fields de- 
stroy superconductivity by polarizing spins, but not by 
inducing vortices. These studies were carried out in a 
very carefully shielded dilution refrigerator to enhance 
the chances that very low heat capacity films would cool. 

The experimental approach used in this work will be 
described in Section II. In Section III the various results 
will be presented. The final Section contains a discus- 
sion of these results along with conclusions that can be 
drawn. 



II. EXPERIMENTAL APPROACH 

These investigations were carried out in a geometry in 
which a SrTi03 (STO) crystal served as both a substrate 
and a gate insulator in a field effect transistor configura- 
tion. To prepare this device, first a small section of the 
unpolished back surface of a 500/im thick single-crystal 
of (100) STO substrate was mechanically thinned 40 ex- 



situ, resulting in this surface and the epi-polished front 
surface being parallel and separated by 45 ± 5 jira. A 
0.5 mm by 0.5 mm, 1000 A thick, Pt "gate" electrode 
was deposited ex situ onto the thinned section of the 
back surface directly opposite the eventual location of 
the measured square of film Platinum electrodes, 100 A 
thick, were also deposited ex situ onto the substrate's epi- 
polished front surface to form a four probe measurement 
geometry. The substrate was then placed in a Kelvinox- 
400 dilution refrigerator/UHV deposition apparatus™ A 
10 A thick under-layer of amorphous Sb and successive 
layers of amorphous Bi (o-Bi) were thermally deposited 
in situ under ultra-high vacuum conditions 10 -9 Torr) 
through shadow masks onto the substrate's front surface. 
The substrate was held at about 7 K during the depo- 
sition process. Films grown in this manner are believed 
to be homogeneously disordered on a microscopic, rather 
than on a mesoscopic scale. 42 

A Keithley 487 voltage source was used to apply volt- 
ages between the film and the gate electrode. In essence, 
the film and the gate electrode formed a parallel plate 
capacitor with the thinned layer of STO serving as the 
dielectric spacer. Applying a positive voltage, Vg, to 
the gate electrode caused electrons to be transferred into 
the film. Since STO crystals have very large dielec- 
tric constants at temperatures below 10 K ( ^ 20,000), 
and since the substrate was greatly reduced in thickness, 
the gate voltage produced large electric fields that fa- 
cilitated large transfers of electrons. Since the dielectric 
constant of STO is known to vary strongly with electric 
field, an analysis was performed that yielded the relation- 
ship between Vq and the areal density of added electrons, 
An, for this substrate and film. At positive gate voltages, 
transferred electron densities were found to be between 
(at V G = 0) and 3.35 x 10 13 cm" 2 (at V G = 42.5 V). 
Application of Vq above 42.5 V did not change any mea- 
surement of film properties relative to those at 42.5 V. 

The sample measurement lines were heavily filtered so 
as to minimize the electromagnetic noise environment of 
the film. The approach was to use RC filters at 300 K 
to attenuate 60 Hz noise, Spectrum Control #1216-001 
7T— section filters at 300 K (in series with 10 kQ resis- 
tors) to attenuate radio frequency noise, and 2 m long 
Thermocoax cables 4 ^ at the mixing chamber stage of the 
refrigerator to attenuate GHz Johnson noise from warmer 
parts of the refrigerator. To avoid complications arising 
from this filtering, measurements were made using DC, 
rather than AC, methods. A 1 nA DC current, I, was 
applied to the films using a Keithley 220 current source. 
Voltage, V, was measured across the centermost 0.5 mm 
by 0.5 mm square of the film using a Keithley 182 volt- 
meter. The sheet resistance of the film, R, was taken to 
be V/I. 

The films clearly fail to cool much below 60 mK, even 
though the dilution refrigerator cools to 7 mK. This is 
almost certainly due to the thermal load on the films 
that is caused by the residual noise environment, together 
with limitations on the thermal grounding of the electri- 
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cal leads. It is very difficult to cool the electrons of a film 
at such low temperatures. Ultrathin films have a negli- 
gible heat capacity and since in this instance, they are 
not immersed in a cooling fluid, the mechanism for the 
electrons to cool is through their thermally (but not elec- 
trically) grounded leads and through contact of the film 
with the "thermal bath" of phonons in the film. At mK 
temperatures, however, electrons and phonons are known 
to decoupled The mixing chamber cools to 7 mK, as 
verified by a 60 Co nuclear orientation thermometer. To 
determine the actual temperature of a film, we used its 
electrical resistance in the insulating regime as a ther- 
mometer. As will be discussed in Section III, a good fit 
to the temperature dependent resistance from 12 K down 
to 60 mK was an exponential activation form character- 
istic of 2D Mott variable range hopping. However, at 60 
mK, R(T) began to deviate from this form and eventu- 
ally became independent of temperature as the mixing 
chamber approached 7 mK. This is consistent with the 
electrons of the film not cooling even while the mixing 
chamber continues to cool. Removal of the 60 Hz filter- 
ing raised the temperature at which the 9.60 A thick 
film began to deviate from the Mott form at 140 mK; 
above this temperature, R(T) was relatively unchanged. 
This suggests that the noise environment prevented the 
film from cooling after the removal of the filter. By an 
extension of this logic, we believe that, with the full fil- 
tering present, the residual noise environment prevented 
the film from cooling below about 60 mK. 

Since in situ rotation of the sample into alignment 
perpendicular to the axis of the magnetic was not pos- 
sible during this series of experiments, Hall effect mea- 
surements that could have yielded the intrinsic charge 
density were not performed. Electrons were stated by 
Buckel to be the charge carriers in a-Bi, from a deter- 
mination of the Hall coefficient in relatively thick (^500 
A ) films. 45 However, the sign of the Hall coefficient is 
known to lead to an erroneous determination of the sign 
of the charge carrier in amorphous semiconductors™ If 
a-Bi acted like a semiconductor in this regard, then holes 
would be the charge carriers, which would make it similar 
to most other superconductors. Thermopower measure- 
ments could accurately ascertain the sign of the of charge 
carriers in a-Bi. 46 Buckel's Hall effect measurements yield 
estimates of areal charge densities of between 10 14 and 
10 15 cm -2 in metallic a-Bi films. 

Careful attention was given to ensure that voltages 
were measured only while the film was in thermal equi- 
librium. Ramping Vq even at slow rates caused heating 
of the film because real currents flowed in the film in re- 
sponse to the displacement currents that were induced by 
changing Vq- Ramping the magnetic field caused Eddy 
current heating. Lastly, it was found that the rate at 
which the film cooled and warmed in high magnetic fields 
decreased dramatically compared to the rate at low fields. 
Thus, the measurement procedure that was used involved 
ramping Vq and/or magnetic field to constant values and 
waiting long times (between three and six hours) to fully 



cool the film to 60 mK. After this, voltages were mea- 
sured as a function of temperature from 60 mK to 1 K 
in controlled steps using temperature feedback and long 
equilibration times (between 5 and 15 minutes per tem- 
perature). Appropriate equilibration times after ramping 
Vet, B, and T were determined by checking reproducibil- 
ity, i.e., by taking measurements after longer wait times, 
with varying rates, and by recording voltage for many 
hours at stable temperatures while watching for drift that 
would signal behavior that was out of equilibrium. 

In a quantum phase transition, the resistances of 2D 
films are expected to obey the finite size scaling func- 
tional form 8 

R/R c = F(\K- K c \/T 1 I v *) (1) 

where F is an unknown function, K is the tuning param- 
eter, K c is the critical value of the tuning parameter, R c 
is the critical resistance, is the correlation length criti- 
cal exponent, and z is the dynamical critical exponent. 
To analyze data using this finite size scaling form, we 
first plotted isotherms of R (K). A crossing point sepa- 
rates the insulating and superconducting phases, yielding 
the critical values R c and K c . As temperature becomes 
too high, subsequent isotherms will not cross at a sin- 
gle point (R C ,K C ), and these temperatures are then ex- 
cluded from the analysis. After determining R c and K C1 
all of the parameters characterizing the scaling analysis 
are known except for the product vz. This product is 
then taken as the unknown, and R/R c is plotted against 
\K -K C \T~ X I VZ for various values of vz. The value of 
vz that produces the best collapse of the data is then the 
exponent product. 



III. RESULTS 

A. The Intrinsic Insulating Regime 

A sequence of a-Bi films was studied. As a func- 
tion of increasing thickness, there were four thicknesses 
that were insulating and nine that were superconduct- 
ing. Here we focus on the four insulating thicknesses 
that we refer to as "intrinsic" insulators. For the data 
shown in Fig. 1, there was no gate bias and no applied 
magnetic field. The temperature dependence of the resis- 
tance, R(T), at each thickness was consistent with Mott 
variable range hopping (VRH) conduction in 2D, as the 
best fits to the data were of the form 

R(T) = R exp([T /T] 1 / 3 ) (2) 

where the constants Ro and To are the resistance pre- 
factor and the activation energy (in K), respectively. 
This is shown in the inset to Fig. 1. For the 9.69 A 
and 9.91 A thick films, data was available to verify this 
form up to 12 K. Unsuccessful attempts were made to 
fit the data with hopping powers of 1/3, 1/2, 0.7 and 1 



3 




0.2 0.4 0.6 0.8 1 
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FIG. 1: R(T) for thicknesses of 9.60, 9.69, 9.91, and 10.22 A 
(top to bottom in the curves in main figure and inset). Data 
for An = 0, B = 0, and T from 60 mK to 1 K, shown as 
R(T) (main) and InR vs.T _1//3 (inset). To and Ro decrease 
with thickness, being 39, 32, 25, and 8 mK and 16406, 16001, 
14317, and 13310 Q , respectively. 
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FIG. 2: i?(T) as a function of An for the 10.22 A thick film 
with B = 0. Data is shown from 60 mK to 1 K. The values 
of An that are shown are 0, 0.62, 1.13, 1.43, 1.61, 1.83, 2.04, 
2.37, 2.63, and 3.35 x 10 13 cm" 2 . Forty four curves of R(T) 
for other values of An are omitted from the plot for clarity. 
Inset: slope of InT from Eq. 3, k vs. An. 



as well as the relationship expected for weak localization 
and electron-electron interactions, 

G = l/R=G + khi(T/To) (3) 

where G is the conductivity and k and To are constants 
Fits were not attempted in which the prefactor R was 
taken to be temperature dependent. 

Mott VRH conduction is found in strongly localized 
systems in the absence of Coulomb interactions between 
electrons. For this work, the absence of Coulomb inter- 
actions may be caused by screening of the electric fields 
in the filmdi because of proximity to the STO substrate, 
which has a very high dielectric constant. 

B. Electrostatic Tuning of the SI Transition in 
Zero Magnetic Field 

In the case of the 10.22 A thick film, which exhibited 
the highest conductivity for the four intrinsically insu- 
lating thicknesses of the sequence, the addition of elec- 
trons to the film induced superconductivity. The evo- 
lution of this electrostatically tuned SI transition with 
charge transfer is shown in Fig. 2. With the charge 
transfer, An = 0, R(T) was well described by Mott hop- 
ping. With An = 0.85 x 10 13 cm -2 , the best fit from 
60 mK to 1 K was that of a ln(T) dependence of the 
conductance on temperature as given in Eq. This was 
the first density at which G(T) was definitely fit better 
by ln(T) than by the Mott hopping form given in Eq. 
[3 At larger values of An, as superconducting fluctu- 
ations became strong at lower temperatures, this ln(T) 
dependence remained the best fit at higher temperatures. 
With An =3.35 x 10 13 cm -2 , the film was fully super- 
conducting (within the scatter of our data) with the su- 
perconducting transition temperature T c = 60 mK. Su- 



perconducting fluctuations were strong between 60 mK 
and about 250 mK, and the conductance was best de- 
scribed by the ln(T) temperature dependence above 250 
mK. The value of the slope of InT is plotted as a function 
of An for superconducting curves in the inset to Fig. 2. 
The fact that the slope varies in a linear fashion with An 
is striking. The actual mechanism for this InT depen- 
dence is unknown. We have previously suggested that it 
may be due to a combination of weak localization and 
electron-electron interactions despite the high value of 
film resistance^ If one assumed that only the electron- 
electron interaction contribution change, then this would 
imply that the Hartree screening parameter changed lin- 
early with An, perhaps contributing to the inducing of 
superconductivity. 

This transition was successfully analyzed using finite 
size scaling employing An as the tuning parameter. This 
suggests that the electrostatically tuned SI transition is 
a quantum phase transition. In the inset to Fig. 3 we 
show R(An) for multiple isotherms between 60 mK and 
140 mK. There is a distinct crossing point at the critical 
electron density, An c = 1.28 x 10 13 cm -2 and the critical 
resistance, R c = 19, 109 Q . In the main part of Fig. 
3, we show the scaling plot. The value of the critical 
exponent product vz that brings about the best collapse 
of the data is 0.7 ± 0.1. The range over which scaling is 
successful is from 60 mK to 140 mK. 

The scaling analysis with uncorrected data appears to 
fail at temperatures above 140 mK because of the InT 
dependence of the conductance in the normal state. It 
is possible to successfully extend the analysis up to 1 
K if this dependence is first removed. This is done by 
assuming that there are two parallel conductance chan- 
nels, one that gives the InT dependence and the other 
that gives superconducting and insulating fluctuations. 
We then find the critical charge transfer to be at An c = 
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FIG. 3: Finite size scaling plot for the 10.22 A thick film with 
B — 0, including data from 60 mK to 140 mK with An as 
the tuning parameter. Fifty four values of An between and 
3.35 x 10 13 cm -2 were included. The best collapse of the 
data was for vz — 0.7 with an uncertainty of ±0.1. Inset: 
R(An) for isotherms between 60 mK and 140 mK. 



FIG. 5: Resistance at 1 K as a function of resistance at 120 
mK for electrostatic (circles) and thickness (squares) tuned 
transitions. As superconductivity develops, the resistance at 
1 K changes much more for thickness tuning than for electro- 
static tuning. 
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FIG. 4: Finite size scaling plot for the 10.22 A thick film with 
B = including values of R(T) from 60 mK to 1 K with An 
as a tuning parameter, having first removed the contributions 
to R(T) from from the conductance channel responsible for 
the InT dependence. The best collapse of the data was for 
vz — 0.7 with an uncertainty of ±0.1. 



0.85 x 10 13 cm -2 since at this temperature R(T) can be 
fit by InT over the entire temperature range (from 60 mK 
to 1 K). The best fit to G(T) at An c is then written as 
GQ C + k c lnT and the critical resistance is the temperature 
independent part of the resistance at An c , or 1/Goc- For 
each An, we subtracted k c lnT from G(T). After this, 
we determined the value of vz that minimized the error 
in the data collapse on the scaling curves in the same 
manner as for the uncorrected data. All values of R(T) 
from 60 mK to 1 K then successfully fell onto the scal- 
ing curves, as shown in Fig. 4. The critical exponent 
product remained at 0.7 ± 0.1. The parallel conductance 
channel responsible for the ln(T) dependence is thus re- 
moved from the analysis using this procedure. This ap- 



proach is more physical than that used by Gantmakher et 
uljlL in the study of perpendicular magnetic field tuned 
SI transitions of In2 03 films, where scaling was success- 
ful at the lowest temperatures, but a positive dR c /dT 
at high temperatures prevented scaling. To broaden the 
temperature range they subtracted a linear temperature 
dependence from R C (T). 

This SI transition appears to involve little change in 
the physical disorder, as determined by noting that the 
resistance at high temperatures (T > IK) changes very 
little as superconductivity develops. This resistance is 
related to the product of the Fermi wave vector, fcp, and 
the electronic mean free path, I. If kpl does not change 
much, then one can assume that disorder does not change 
much. In Fig. 5, we show the resistance at low temper- 
ature (120 mK) as a function of the resistance at 1 K 
for the electrostatic tuned transition and a previously 
reported thickness-tuned transition. 25 One can see that, 
as insulating behavior changes to superconducting be- 
havior, the thickness-tuned transition takes place with a 
much larger change in the high temperature resistance 
than does the electrostatically tuned system. To fur- 
ther the idea that electrostatic doping adjusts electronic 
properties but not disorder, we note that for intrinsically 
superconducting films with transition temperatures be- 
tween 1 and 4 K, the application of Vq = ±50 V shifted 
transition temperatures by as much as 40 mK, but these 
shifts were accompanied by very small changes in the re- 
sistance at 10 K. The later changed by no more than 
10 — 20 Q with normal state resistances of 3-5 kQ. 



C. Electrostatic Tuning of the SI Transition in a 
2.5 T Parallel Magnetic Field 

Electrostatic tuning of the SI transition of the 10.22 A 
thick film was also carried out in a parallel magnetic field 
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FIG. 6: R{T) as a function of An for the 10.22 A thick film 
with B = 2.5 T. Values of An shown are 0, 0.74, 1.28, 1.57, 
1.87, 2.30, 2.74 and 3.13 x 10 13 cm" 2 , from top to bottom. Fif- 
teen curves of R(T) at other values of An have been omitted 
from the plot for clarity. 



FIG. 7: Finite size scaling plot for the 10.22 A thick film with 
B = 2.5 T including values of R(T) from 100 mK to 200 mK 
with An as the tuning parameter. Twenty three values of 
An between and 3.13 x 10 13 cm 
collapse of the data was for vz - 
of ±0.1. Inset: R(An) for isotherms between 100 mK and 
200 mK, exhibiting a well-defined crossing point at An = 1.7 
xl0 13 cm- 2 and R = 18,300 Q. 



2 were included. The best 
0.65 with an uncertainty 



of 2.5 T. The data is shown in Fig. 6. The application of 
magnetic field to an ungated insulating film increased its 
activation energy. Adding electrons increased the film's 
conductivity and, at high enough density, induced a tran- 
sition to superconductivity. 

Comparing Figs. 2 and 6, a major difference between 
the electrostatically tuned SI transitions in zero and non- 
zero magnetic fields is evident: at low temperatures, in 
nonzero field, a regime exists in which the resistance 
at the lowest temperatures is larger than that expected 
from the extrapolation of R(T) from higher tempera- 
tures. This is most obvious for electron densities around 
the critical density, where curves that appear to be head- 
ing towards zero resistance as temperature is lowered sud- 
denly undergo a change in slope and appear to be insulat- 
ing when extrapolated to the limit of zero temperature. 

A finite size scaling analysis was carried out for the 
electrostatically tuned transition in field. In the inset to 
Fig. 7 we show R(An) for isotherms between 100 and 
200 mK. A distinct crossing point is apparent, yielding 
An c = 1.7 x 10 13 cm -2 and R c = 18,300 tt. In the main 
body of Fig. 7, we show the scaling plot. The value of 
the exponent product that minimized the collapse of the 
data was 0.65 ± 0.1, which is slightly lower than the elec- 
trostatic tuned transition in zero field but is in agreement 
with it if the uncertainty of the analysis is taken into ac- 
count. Data below 100 mK does not fall onto the scaling 
plot because of the above mentioned excess resistance. 

The successful scaling analysis implies that the sys- 
tem is moving towards insulating and superconducting 
ground states that are separated by a quantum critical 
point. The break down in scaling implies that instead of 
reaching this insulating ground state, one with a higher 
resistance is actually achieved. This will be discussed in 
more detail in in Section F. 
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FIG. 8: R(T) as a function of B for an insulating film with 
superconductivity induced by a charge transfer An = 3.35 x 
10 13 cra~ 2 . Values of B are 2 (bottom), 2.5, 3.5, 4.25, 4.375, 
4.75, 5, 5.75, 6.5, 8, 9, and 11 T (top). Seven R(T) curves for 
other values of B have been omitted from the plot for clarity. 



D. Parallel Magnetic Field Tuning of the SI Tran- 
sition 

In the case of the 10.22 A thick film, after supercon- 
ductivity was induced by adding an areal density of car- 
riers An = 3.35 x 10 13 cm -2 , this induced supercon- 
ductivity was quenched by a parallel magnetic field, and 
a detailed study of the parallel magnetic field tuned SI 
transition was conducted. In Fig. 8, we show this transi- 
tion. With B = 0, the superconducting transition tem- 
perature, T c = 60 mK (taken to be the highest temper- 
ature at which resistance is zero within the scatter due 
to noise), the mean field transition temperature (taken 
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FIG. 9: Finite size scaling plot for the 10.22 A thick film with 
An = 3.35x 10 13 cm -2 with B as the tuning parameter for 
150 mK < T < 340 mK. The best collapse of the data was 
for vz — 0.65 with an uncertainty of ±0.1. Inset: Isotherms 
of R(B) at temperatures between 150 mK and 340 mK. 



FIG. 10: Phase diagram of B c vs. (An - An c ) for the 10.22 
A thick film. The best fit is by a power law with an exponent 
of 0.33. 



to be the temperature at which resistance is half of the 
normal state peak) is 90 mK, and the peak in R(T) be- 
low which there are strong superconducting fluctuations 
is at 250 mK. The normal state, taken to be at tempera- 
tures in excess of that at which R(T) exhibits a peak, is 
best described by a conductance that has a logarithmic 
temperature dependence. This superconductor at B = 
becomes an insulator at high fields that is again well de- 
scribed by 2D Mott VRH at temperatures greater than 
130 mK. 

Analysis of this transition using finite size scaling is 
successful over a wide range of temperature, provided 
that the regime of excess resistance discussed above is 
excluded from the analysis. The scaling analysis is shown 
in Fig. 9. The inset to Fig. 9 shows a set of isotherms 
of resistance vs. field between 150 and 340 mK, which 
exhibits a distinct crossing point at R c = 17, 285 and 
B c = 4.625 T. The value of the exponent product vz 
that minimized collapse of the data is 0.65 ± 0.1. Within 
the experimental uncertainty, this value is the same as 
that found for the electrostatically tuned transitions in 
zero field and 2.5 T. 



E. Phase Diagrams 

We have carried out a systematic study of the parallel 
magnetic field tuned SI transition at various strengths of 
electrostatically induced superconductivity. The scaling 
procedure for the B-tuned transition for An = 3.35 x 10 13 
cm -2 discussed in the previous section was repeated at 
An =1.66, 2.25, and 2.80 x 10 13 cm -2 , yielding the same 
exponent product vz = 0.65 ±0.1. 

At values of An = 1.83, 2.00, 2.40, 2.59 and 2.96 x 
10 13 cm -2 , R(B) was measured at 100 and 120 mK and 
the crossing points determined. These temperatures were 
within the range used to determine the crossing points 



when a full scaling analysis was carried for An = 1.66, 
2.25, and 2.80 x 10 13 cm -2 using data at many different 
temperatures. Thus, they are greater than the temper- 
atures at which excess resistance occurs. Because the 
crossing points were well defined when there was exten- 
sive data, one can be confident that the intersection at 
two temperatures would be a reliable determination of 
the critical resistance and critical field. 

Using both the full scaling analyses and the reduced 
crossing-point analyses, we were able to map out the 
variations of both B c and R c with An as shown in Figs. 
10 and 11, respectively. The critical field, B c , increased 
with An, as (An — An c ) 0,33 . Qualitatively, this is what 
would be expected, since adding carriers strengthens su- 
perconductivity and destroying stronger superconductors 
would require higher magnetic fields. The critical resis- 
tance decreased linearly with An. The data from the 
electrostatically tuned transition in a magnetic field of 
2.5 T is included on these plots. If this data is excluded, 
the functional forms are unchanged. 

The exponent product vz is the same, within experi- 
mental uncertainty, for the parallel magnetic field tuned 
transitions at various strengths of electrostatically in- 
duced superconductivity and the electrostatically-tuned 
transitions in zero and finite magnetic fields. This sug- 
gests that the quantum phase transitions belong to the 
same universality class and that these phase diagrams de- 
mark the insulating and superconducting regimes with a 
robust phase transition line that can be crossed in either 
direction, by tuning B or An. This result is differs from 
that found by Markovic et alM in which the exponent 
products of thickness and perpendicular field tuned tran- 
sitions were found to be different. Note that the regime 
of excess resistance found in parallel magnetic fields may 
imply that the insulating ground state may be different 
from that which might be deduced from the analysis of 
data obtained at high temperatures and used to construct 
Figs. 11 and 12. 
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0.5 1 1.5 2 

An-An c (10 13 cm" 2 ) 

FIG. 11: Phase diagram of R c vs. (An - An c ) for the 10.22 
A thick film. The best description is a linear relation. 



F. Excess Resistance in High Magnetic Fields 

In films whose parameters place them on the bor- 
der between the insulating and superconducting regimes, 
or well into the insulating regime, we observed behav- 
ior in which resistance in magnetic fields became higher 
at low temperatures than what might have been ex- 
pected from extrapolation of R(T) at higher tempera- 
tures. This regime appears to involve physics that is 
different from that which determines quantum critical- 
ly. While a scaling analysis is successful including data 
obtained at higher temperatures, it fails when data from 
this regime is included. 

To illustrate this excess resistance, we now discuss in 
more detail the data on the parallel magnetic field-tuned 
transition shown in Fig. 9. Here, the resistance is larger 
than what is expected by extrapolating the curves of 
R(T) down from temperatures in excess of 150 mK. Be- 
tween 3 T and the parallel critical field of 4.625 T, while 
superconducting fluctuations cause R(T) curves to head 
towards zero resistance as temperature is lowered, there 
is a small upturn in resistance as 60 mK is approached. 
This upturn in R(T) is found at all higher fields in the 
insulating regime. It is easily visible at fields between 
4.625 T and 7 T, where the upturn occurs simultane- 
ously with a minimum in R(T). In fields in excess of 7 
T, there are no minima in R(T) since the general shape 
is that of an insulating curve with dR/dT < for all T, 
but there is still an excess resistance at low temperatures. 
Here, resistance is in excess of an extrapolation of Mott 
VRH, which describes the data from 130 mK up to 1 K. 
To illustrate this, R(T) in fields of 8, 9, and 11 T from 
Fig. 8 are re-plotted in Fig. 12 in the form ln(R) vs. 
T -1 / 3 . At higher temperatures, 2D Mott VRH is seen to 
be a good fit, while below about 130 mK, the resistance 
becomes greater than extrapolations using the Mott form 
from higher temperatures. 

The minima in R(T) and the increase in resistance 
at temperatures below those at which the minima occur 



10.2 rr 




T -w (K -ifl) 

FIG. 12: With superconductivity first induced in the 10.22 
A thick film with An = 3.35 x 10 13 cm" 2 , ln(R) vs. T" 1/3 
for insulating curves induced by 8 (bottom), 9, and 11 (top) 
T. The data is fit well by Mott VRH for T >150 mK, but 
deviates to a higher resistance below 150 mK. Straight lines 
have been added as guides to the eye. 



are both qualitatively similar to features in granular su- 
perconducting films. In granular films, minima in R(T) 
and the subsequent increase of resistance with decreasing 
temperature are believed to be indicative of persistence 
of superconductivity on mesoscopic sized grains, where 
inter-granular transport is via single-particle tunneling. 
We have not observed reentrance in the absence of mag- 
netic field, either in the electrostatically tuned SI transi- 
tion or the thickness-tuned transition (see Figs. 1 and 2). 
This suggests that there are no mesoscopic scale clusters 
in the film. 

The magnitude of the excess resistance in magnetic 
field in these films is much smaller than that in granular 
films. For instance, for B = 8, 9, and 11 T, resistances 
become about 4-8 % higher at 65 mK than one would 
expect from extrapolation. At these fields, the minima 
occur at about 140 mK. In contrast, in granular films at 
temperatures a factor of two lower than the temperature 
at a minimum, resistance is found to be some 16 times 
higher than is predicted, due to a (T/T m ^ n ) 4 dependence 
associated with the opening of the energy gap as super- 
conductivity develops on the grains 

The highest temperature at which excess resistance is 
observed in a given field, denoted as T onsetl increases 
with increasing field. This was checked for the four sets 
of parallel field induced SI transitions for which full scal- 
ing analyses were performed, as well as for an intrinsically 
insulating film where superconductivity was not induced 
electrostatically. For each superconducting film, the on- 
set of excess resistance appears around the critical mag- 
netic field of the SI transition. As mentioned before, B c 
increases with An. Interestingly, these T onset (B) curves 
collapse if shifted horizontally to align the critical fields. 
This is shown in Fig. 13. This implies that when super- 
conductivity is stronger (at larger An) it takes a larger 
field to induce excess resistance than it takes when su- 
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FIG. 13: T onset vs. B — B c for superconducting films with 
various values of An. An = 1.66 (diamonds), 2.25 (triangles), 
2.80 (squares), and 3.35 (circles) x 10 13 cm -2 . The data has 
been collapsed by shifting these curves by 2, 2.7, 3.5, and 3.75 
T, respectively, from the unadjusted values. 



perconductivity is weaker (at smaller An) . It also implies 
that once excess resistance has been induced, it develops 
with increasing field in a manner that is independent of 
the initial strength of superconductivity. 

There have been several observations of a peak in 
the magnetoresistance at fields in excess of the criti- 
cal field for the SI transition in thicker films (^ 50 — 
300 A) of ln 2 3 and of TiN, both in fields aligned 
perpendicula u 13 i 30 i 31 i 33 i 34 i 35 to and parallel 38 to the film 
plane. This peak has even been observed in insulating 
Be films Typically, this resistance peak is found to be 
larger in films with high normal state resistances (and 
correspondingly low values of T c ) than in those with low 
normal state resistances and higher values of T c & When 
we induce superconductivity in the 10.22 A thick film 
with An = 3.35 x 10 13 cm. -2 , with R N - 16 k^ and 
T c ^ 60 mK, we do not observe a peak in resistance as 
a function of field. In Fig. 14, we show R(B) for various 
isotherms between 65 mK and 200 mK. For all isotherms, 
dR/dB > at high fields. This appears to be a general 
feature of our films, as we observed dR/dB > up to 12 
T at all temperatures in the other parallel magnetic field 
tuned transitions at smaller values of An. 



FIG. 14: R(B) for the 10.22 A thick film with An = 3.35 x 
10 13 cm -2 . The isotherms correspond to 65, 75, 85, 100, 120, 
150, and 200 mK. At high fields, dR/dB > 0. 



deposited on top of 10 A thick a-Sb layer on a STO sub- 
strate that had not been thinned. This film was super- 
conducting with a transition temperature T c = 130 mK 
and a normal state sheet resistance, Rn ^ 9500 Q. Paral- 
lel magnetic fields induced insulating behavior and finite 
size scaling was successful over a range of temperatures 
from 150 to 400 mK. A distinct crossing point was found, 
with R c = 10,460 n and B c = 9.18 T. The critical ex- 
ponent product that produced the best collapse of the 
data was 0.75 ± 0.1, slightly higher than found for the 
field-tuned transition of an electrostatically induced su- 
perconducting film, but agreeing within the uncertainty. 
Excess resistance was observed in fields ranging from 8 T 
to 12 T at temperatures below about 150 mK. At 12 T, 
a fit to R(T) by Mott VRH was successful down to 150 
mK. For data below 150 mK, the excess resistance again 
prevented successful scaling. For T > 400 mK, the neg- 
ative slope of dR c /dT prevented successful scaling. At 
high fields and at all temperatures, dR/dB was always 
positive. 



IV. DISCUSSION 

A. Magnetic Field Alignment 



G. Comparison of the Field- Tuned SI Transitions 
of Intrinsic and Electric-Field Induced Superconduc- 
tors 

A separate experiment on parallel field tuning of the 
SI transition was conducted on an intrinsically supercon- 
ducting film. The results of a finite size scaling analysis 
were nearly identical to those found for the transition of 
an electrostatically-induced superconducting film. The 
excess resistance in high fields was also nearly identical 
to that found in an electrostatically-induced supercon- 
ducting film. The sample was a 10.25 A thick a-Bi film 



An important technical issue for studies in parallel 
magnetic fields is the alignment of the plane of the sub- 
strate with the field. Our estimate of the maximum angu- 
lar misalignment of the is 1°. This is based on multiple 
geometrical constraints. For the magnetic field tuned 
SI transition for superconductivity induced by adding 
An = 3.35 x 10 13 cm -2 , the transition temperature is 60 
mK, the "mean field" transition temperature is 90 mK, 
and the critical field is 4.625 T. This critical field is 27 
times higher than the limiting value necessary to destroy 
superconductivity by aligning spins, found by Clogston 
and Chandrasekhar to be 1.9 T/K. An enhancement of 
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this magnitude might be expected as Bi is a heavy metal 
with strong spin-orbit interaction and consequently a-Bi 
would be expected to have a short spin orbit scattering 
time. At each electron scattering event, the spin flips, 
leading to a substantially enhanced value of the criti- 
cal field relative to the Clogston/Chandrasekhar limiting 
field. If we were to assume that this result was due to 
a perpendicular field component resulting from misalign- 
ment larger than our estimate, a misalignment of about 
7 to 11 degrees would be necessary to produce the crit- 
ical fields in Fig. 10. We base this on comparison with 
work on similar a-Bi films by Markovic et alML in which a 
perpendicular field of 0.6 T was needed to quench super- 
conductivity in a film that has a transition temperature 
between that found here at An = 1 x 10 13 cm -2 and 3.35 
x 10 13 cm -2 , based on its shape. This is unreasonably 
large, given the geometrical constraints. 



B. SI Transitions 

In aggregate, finite size scaling analyses of the vari- 
ous SI transitions tuned by thickness, electrostatic elec- 
tron doping, and perpendicular and parallel magnetic 
fields have yielded a wide range of results. Scaling 
analyses of measurements in which perpendicular mag- 
netic field^i 18 ! 19 ! 31 ! 38 ! 50 or thickness 20 have been tuned 
have usually resulted in vz ^ 1.3 and z = 1, which is 
consistent both with the scaling theory, 9 (2+1) dimen- 
sional XY model with disorder^ and other and numeri- 
cal studiesj 52 i 53 i 54 i 55 i 56 i 57 This value has been suggested 
to be consistent with percolation*22i2£ An exception is 
the case of a-Bi, in which tuning the transition by a per- 
pendicular field has yielded z ^ 0.7, while tuning with 
thickness has yielded z ^ 1.3 Assuming that z = 1, this 
value disagrees with what was believed to be an theorem 
that predicts v > 1 in two dimensions in the presence 
of disorder. 58 The possible irrelevance of this theorem 
to the SI transition has been discussed most recently by 
Chamon and Nayak™ The present measurements sug- 
gest that transitions tuned by electron doping and par- 
allel magnetic fields have exponent products of 0.65 to 
0.75. Thus if z = 1, vz = 0.7 is consistent with earlier 
perpendicular field scaling, as well as numerical studies 
of the 3D XY model 7-51 and the Boson-Hubbard mod- 
els in the absence of disorder™ This is suggestive that 
thickness-tuned transitions in a-Bi may be percolative, 
while transitions for other tuning parameters are not. 

The values of critical resistances for all of the transi- 
tions of electrostatically induced superconductivity are a 
factor or two or three higher than those found in other SI 
transitions and the universal value of 6455 Vl predicted by 
the dirty Boson model. While the critical resistance has 
been found experimentally to be no n- universal ^LA& it has 
been studied in various materials and samples with dif- 
ferent levels of disorder. Here, in a sample with a static 
level of disorder, we have shown the critical resistance 
decreases linearly with increasing electron density. 



Excess resistance appears to be a low temperature fea- 
ture of films in the insulating regime in the presence of 
parallel magnetic fields. Fields have produced this both 
insulating regime when superconductivity was quenched 
by field and in the intrinsic insulating state. The features 
of R(T) in parallel magnetic fields are qualitatively simi- 
lar to features in R(T) in zero field for granular films. The 
conventional view is that quench-condensed films grown 
with either an a-Ge or an a-Sb underlayers are disordered 
on an atomic rather than a mesoscopic scale. Since these 
films anneal at temperatures around 25 K, whereupon 
they change structurally into semimetals, we cannot ex- 
amine them structurally to determine the level of homo- 
geneity. However, since we have not observed reentrance 
in the absence of magnetic field, either in the electrostat- 
ically tuned SI transition or the thickness-tuned transi- 
tion, we believe that there are no mesoscopic scale clus- 
ters in the film. It is possible that the parallel magnetic 
field induces spatial inhomogeneity of the amplitude of 
the superconducting order parameter that mimics clus- 
tering. 

In the present work, we found that the relation be- 
tween critical parallel magnetic field and electron den- 
sity is B c ^ (An - An c ) 33 , whereas Markovic et alM 
found the relationship between perpendicular magnetic 
field and film thickness, d, as B c ^ (d — d c ) 0,7 for perpen- 
dicular fields. There appears to be no theory constraining 
the relationship between An and B c . It appears that the 
line of criticality can be crossed changing either B or An 
with a finite size scaling analysis yielding the same ex- 
ponent product, whereas in the work of Markovic et al. 
tuning with d and with B produced different exponent 
products. 

Flattening of R (T) below about 60 mK in zero mag- 
netic field may only be a consequence of failure to cool the 
film. In nonzero magnetic field its occurrence at tem- 
peratures higher than 60 mK may be due to enhanced 
Eddy current heating or there may be an intrinsic metal- 
lic regime. 

Interestingly, we have not observed a peak in the mag- 
netoresistance after superconductivity has been quenched 
by magnetic fields, even though our transition tempera- 
tures are very low and normal state resistances are very 
high, which is the regime in which the largest peaks of 
In203 films were found™ We do not find a peak ei- 
ther below 150 mK, in the regime of excess resistance, 
or above this temperature, where Mott VRH fits well. 
Though our fields are aligned parallel to the plane of 
the film, large peaks in magnetoresistance were found by 
Gantmakher et al. in alignments both perpendicular and 
parallel to the film™^ Though there is theoretical work 
suggesting that this metallic regime is intrinsic, 26 there 
have been suggestions that this is a consequence of in- 
homogeneity in the films We suggest that further 
experimental studies of the homogeneity of films exhibit- 
ing peaks in magnetoresistance in the presence of large 
magnetic fields would help resolve the issue of whether 
these magnetoresistance peaks are intrinsic properties of 
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structurally and chemically homogeneous films. 

V. CONCLUSIONS 

We have investigated the two-dimensional 
superconductor-insulator transition in disordered 
ultrathin a-Bi films by use of electrostatic electron 
doping using the electric field effect and by the use of 
parallel magnetic fields. Electrostatic doping was carried 
out in both zero and nonzero magnetic fields, and 
magnetic tuning was conducted at multiple strengths of 
electrostatically induced superconductivity. The various 
transitions were analyzed using finite size scaling to 
determine the critical exponent products of the quantum 
phase transitions, which were all found to be vz = 0.65 
or 0.7 ±0.1. The critical parallel magnetic field increased 
with electron transfer, An, as (An c — An) 0,33 , while 
the critical resistance decreased linearly with electron 
transfer. 

An anomalous regime of excess resistance that is in- 
duced by a parallel magnetic field was also observed. This 
excess resistance was not observed in zero magnetic field 
for either thickness or electrostatic tuning, the existence 
of this regime necessitated that data below 100 to 150 
mK be excluded in order for scaling to be successful. 

Although there is a long history of experimental and 
theoretical investigation of the control of superconduc- 
tivity with electric fields beginning with Glover and 
Sherill, 61 there has been no significant tuning of super- 
conductivity in metallic systems, and models based on 
the BCS theoryfiSifiiSi only treat the effect on the tran- 
sition temperature of changes in the density of states in 
response to changes in the carrier concentration. They 
do not include an issue that is relevant here, the appar- 
ent insulator to metal transition that accompanies the 
superconductor to insulator transition^ The change in 
the coefficient of the InT behavior of the resistance in the 
normal state reported in Parendo et al 39 implies that the 
effective electron-electron interaction is also changed as 
the carrier concentration (electron density) is increased. 
This is not included in any of the theoretical treatments 
and is a challenge for theory. 

The parallel field tuned SI transition leads to values of 
the critical field that are the order of a factor of 30 in ex- 
cess of the Pauli limiting field for a superconductor with 
a transition temperature of 60 mK. This is to be expected 



for a system with strong spin-orbit scattering. In this 
instance the rather extensive theoretical treatments^ of 
the critical fields of superconductors are not applicable as 
quantum fluctuations and the fact that the normal phase 
is an insulator and not a metal have not been taken into 
account. 

Finally the excess resistance in the insulating phase of 
the parallel field tuned transition, which actually appears 
like quasi-rentrant superconductivity, strongly suggests 
the existence of a regime in which the superconducting 
order parameter is inhomogeneous, i.e., that there may 
be superconducting droplets in an insulating matrix. Al- 
though this was observed only in the field-tuned tran- 
sition and not in the charge-tuned transition, the idea 
is fairly generic in theories of SI transitions in homoge- 
neously disordered films. The case of perpendicular fields 
has been considered by Spivak and Zhou§& and by Gal- 
itski and Larkin 12 . For thin films with strong spin-orbit 
scattering, as would be the case for a-Bi, Zhou££ has 
shown that there should be a glassy phase. We have ob- 
served very slow relaxation in hgih magnetic fields, which 
could be evidence of such a reigme. In all instances 
the inhomogeneous regime is confined to a narrow re- 
gion around the critical field. Recently Skvortsov and 
Feigelman— have considered the suppression of super- 
conductivity by disorder in films of high dimensionless 
conductance. They find an inhomogenous regime close 
to the critical conductance where superconductivity is 
completely suppressed. This is not seen in our electro- 
statically tuned transition, perhaps due to the relatively 
low normal state conductance in our films. 
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